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Abstract 

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sen- 
r^. sitivity to the physics of jet quenching. In a sample of lead-lead collisions at ^/snn = 2.76 TeV 
corresponding to an integrated luminosity of approximately 7 fib' 1 , ATLAS has measured jets with 
a calorimeter over the pseudorapidity interval \r]\ < 2.1 and over the transverse momentum range 
38 < p T < 210 GeV. Jets were reconstructed using the anti-/c t algorithm with values for the distance 
parameter that determines the nominal jet radius of R = 0.2,0.3,0.4 and 0.5. The centrality depen- 
dence of the jet yield is characterized by the jet "central-to-peripheral ratio," R C p- Jet production is 
found to be suppressed by approximately a factor of two in the 10% most central collisions relative 
to peripheral collisions. R C p varies smoothly with centrality as characterized by the number of partic- 
ipating nucleons. The observed suppression is only weakly dependent on jet radius and transverse 
momentum. These results provide the first direct measurement of inclusive jet suppression in heavy 
ion collisions and complement previous measurements of dijet transverse energy imbalance at the 
LHC. 
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Abstract 

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to 
the physics of jet quenching. In a sample of lead-lead collisions at -^/snn = 2.76 TeV corresponding to 
an integrated luminosity of approximately 7 fjh , ATLAS has measured jets with a calorimeter over the 
pseudorapidity interval \rj\ < 2.1 and over the transverse momentum range 38 < pt < 210 GeV. Jets 
were reconstructed using the anti-fc t algorithm with values for the distance parameter that determines the 
nominal jet radius of R = 0.2,0.3,0.4 and 0.5. The centrality dependence of the jet yield is characterized 
by the jet "central-to-peripheral ratio," Rqp- Jet production is found to be suppressed by approximately a 
factor of two in the 10% most central collisions relative to peripheral collisions. Rcp varies smoothly with 
centrality as characterized by the number of participating nucleons. The observed suppression is only weakly 
dependent on jet radius and transverse momentum. These results provide the first direct measurement of 
inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse 
energy imbalance at the LHC. 
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1. Introduction 

Collisions of lead ions at the LHC are expected 
to create strongly interacting matter at the high- 
est temperatures ever produced in the laboratory 
1\. This matter may be deconfined with a high 
density of unscreened colour charges. High trans- 
verse momentum (pt) quarks and gluons generated 
by hard-scattering processes have long been consid- 
ered an important tool for probing the properties of 
the matter created in ultra-relativistic nuclear col- 
lisions. The energy loss of the partons propagating 
through the matter may provide direct sensitivity 
to the colour charge density and to the transport 
properties of the matter (2H1] . Indirect observations 
of substantial parton energy loss or "jet quench- 
ing" via suppressed single high-px hadron yields 
[SH5] and disappearance of the dijet contribution 
to di-hadron correlations [TU] have contributed 
to the conclusion that Au+Au collisions at RHIC 
produce a quark-gluon plasma [11] [12]. Observa- 
tions of highly asymmetric dijets in central Pb+Pb 
collisions at the LHC [T3TH5] can be understood in 



the context of "differential" jet quenching, where 
one parton produced from an initial hard-scattering 
loses significantly more energy than the other, pos- 
sibly as a result of different path lengths of the par- 
tons in the matter 16J. However, the asymmetry 
is not sensitive to situations where the two jets in 
a dijet pair lose comparable amounts of energy, so 
other measurements are required to probe "inclu- 
sive" jet quenching. 

The inclusive, per-event jet production rate pro- 
vides such a measurement. Energy loss of the par- 
ent partons in the created matter may reduce or 
"suppress" the rate for producing jets at a given 
transverse momentum. Such energy loss is expected 
to increase with medium temperature and with in- 
creasing path length of the parton in the medium 
|17j . As a result, there should be more suppression 
in central Pb+Pb collisions, which have nearly com- 
plete overlap between the incident nuclei, and lit- 
tle or no suppression in peripheral collisions where 
the nuclei barely overlap. In the absence of energy 
loss, the jet production rate is expected to vary 
with Pb+Pb collision centrality approximately in 
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proportion to iVcoli> the number of nucleon-nucleon 
collisions that take place during a single Pb+Pb 
collision. The jet suppression may be quantified 
using the central-to-peripheral ratio, Rcp, the ra- 
tio of the per-event jet yields divided by the number 
of nucleon-nucleon collisions in a given centrality 
bin to the same quantity in a peripheral centrality 
bin. The quantity, Rcp > has the advantage that po- 
tentially large systematic uncertainties, especially 
those arising from systematic errors on the jet en- 
ergy scale, largely cancel when evaluating the ratios 
of jet spectra within the same data set. The varia- 
tion of the suppression with jet transverse momen- 
tum and with collision centrality will depend both 
on the energy loss mechanism and on the experi- 
mental definition of the jet. In the case of radiative 
energy loss, jet energies can be reduced by greater 
"out-of-cone" radiation, which should be more se- 
vere for smaller jet radii (T8 20J . Thus, a measure- 
ment of Rqp for different jet radii could clarify the 
role of radiative energy loss in jet quenching and 
probe the angular distribution of medium-induced 
radiation. 

This Letter presents measurements of the inclu- 
sive jet Rqp in Pb+Pb collisions at a nucleon- 
nucleon centre-of-mass energy of = 2.76 TeV 
using data collected during 2010 corresponding to 
an integrated luminosity of approximately 7 fjh . 
Results are presented for jets reconstructed from 
energy deposits measured in the ATLAS calorime- 
ters using the anti-fc t jet-finding algorithm 2 Ij . The 
anti-fc t reconstruction was performed separately for 
four different values of the anti-fct distance param- 
eter, i?, that specifies the nominal radius of the re- 
constructed jets, R = 0.2,0.3,0.4 and 0.5. For the 
remainder of the Letter the term "radius" will refer 
to the distance parameter, R. The assumed addi- 
tive contribution of the underlying event to each 
jet was subtracted on a per-jet basis, and the i?cp 
values were calculated after unfolding the jet spec- 
tra for distortions due to intrinsic jet resolution and 
underlying event fluctuations. 

2. Experimental setup and trigger 

The measurements presented here were per- 
formed using the ATLAS calorimeter, inner de- 
tector, trigger, and data acquisition systems |22| . 
The ATLAS calorimeter system consists of a liq- 
uid argon (LAr) electromagnetic (EM) calorime- 
ter covering I77I < 3.2, a steel-scintillator sam- 
pling hadronic calorimeter covering < 1.7, a 



LAr hadronic calorimeter covering 1.5 < \rj\ < 3.2, 
and two LAr electromagnetic and hadronic forward 
calorimeters (FCal) covering 3.2 < |?7| < 4.£H. The 
hadronic calorimeter granularities or cell sizes in 
Ari x are 0.1 x 0.1 for |r?| < 2.5 and 0.2 x 0.2 for 
2.5 < |?7| < 4.9^ The EM calorimeters are longitu- 
dinally segmented into three compartments with an 
additional pre-sampler layer. The EM calorimeter 
has a granularity that varies with layer and pseu- 
dorapidity, but which is generally much finer than 
that of the hadronic calorimeter. The middle sam- 
pling layer, which typically has the largest energy 
deposit in EM showers, has a Ai] x A<p granularity 
of 0.025 x 0.025 over \r)\ < 2.5. 

Charged particles associated with the calorime- 
ter jets were measured over the pseudorapidity in- 
terval |f?| < 2.5 using the inner detector [23 . The 
inner detector is composed of silicon pixel detec- 
tors in the innermost layers, followed by silicon mi- 
crostrip detectors and a straw-tube tracker, all im- 
mersed in a 2 T axial magnetic field provided by 
a solenoid. Minimum bias Pb+Pb collisions were 
identified using measurements from the zero-degree 
calorimeters (ZDCs) and the minimum-bias trig- 
ger scintillator (MBTS) counters. The ZDCs are 
located symmetrically at z = ±140 m and cover 
|?7| > 8.3. In Pb+Pb collisions the ZDCs primar- 
ily measure "spectator" neutrons - neutrons from 
the incident nuclei that do not interact hadroni- 
cally. The MBTS measures charged particles over 
2.1 < \r\\ < 3.9 using two sets of counters placed at 
z = ±3.6 m. Events used in this analysis were se- 
lected for recording by the data acquisition system 
using a logical or of ZDC and MBTS coincidence 
triggers. The MBTS coincidence required at least 
one hit in each side of the detector, and the ZDC co- 
incidence trigger required the summed pulse height 
from each calorimeter to be above a threshold set 
below the single neutron peak. 



1 ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the y 
axis points upward. Cylindrical coordinates (r, <j>) are used 
in the transverse plane, <f> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of 
the polar angle 9 as r) = — lntan(6>/2). 

2 An exception is the third (outermost) sampling layer, 
which has a segmentation of 0.2 X 0.1 up to \r)\ = 1.7. 
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3. Event selection and centrality definition 

In the offline analysis, Pb+Pb collisions were re- 
quired to have a primary vertex reconstructed from 
charged particle tracks with pr > 500 MeV. The 
tracks were reconstructed from hits in the inner 
detector using the standard ATLAS track recon- 
struction algorithm [53] with settings optimized for 
the high hit density in heavy ion collisions [25] . 
Additional requirements of a ZDC coincidence, at 
least one hit in each MBTS counter, and a differ- 
ence in times measured by the two sides of the 
MBTS detector of less than 3 ns were imposed. 
The combination of the ZDC and MBTS condi- 
tions and the primary vertex requirement efficiently 
eliminates both beam-gas interactions and photo- 
nuclear events [25]. These event selections yielded 
a total of 51 million minimum-bias Pb+Pb events. 
Previous studies [35] indicate that the combination 
of trigger and offline requirements select minimum- 
bias hadronic Pb+Pb collisions with an efficiency 
of 98 ± 2%. 

The centrality of Pb+Pb collisions was charac- 
terized by ££^ Cal , the total transverse energy mea- 
sured in the forward calorimeters. The distribution 
of T,E^ Cal was divided into intervals correspond- 
ing to successive 10% percentiles of the full cen- 
trality distribution after accounting for the missing 
2% most peripheral events. A standard Glauber 
Monte-Carlo analysis [371 HE] was used to esti- 
mate the average number of participating nucleons, 
(iVpart), and the average number of nucleon-nucleon 
collisions, (N co n), for Pb+Pb collisions in each of 
the centrality bins. The results are shown in Ta- 
ble [I] The Rqp measurements presented here use 
the 60-80% centrality bin as a common peripheral 
reference. The Rcp calculation requires the ratio, 



Table 1: Results of Glauber model evaluation of (iVpart) 
and associated errors, (iV ca ii) , the N co u ratios, -R co ii, and 
fractional errors on R co \\ for the centrality bins included in 
this analysis. 



Centrality 


(Apart) 


(JVcou) 


-Rcoll 


- 10% 


356 + 2 


1500 + 115 


57+11% 


10 - 20% 


261 + 4 


923 ± 68 


35 ± 11% 


20 - 30% 


186 + 4 


559 ± 41 


21 ± 9% 


30 - 40% 


129 + 4 


322 ± 24 


12 ± 8% 


40 - 50% 


86 + 4 


173 ± 14 


6.5 ± 6% 


50 - 60% 


53 + 3 


85 + 8 


3.2 ± 4% 


60 - 80% 


23 + 2 


27 + 4 


1 



R coll = (N con )/(N™-W), where (A C 6 V°) is the av- 
erage number of collisions in the 60-80% centrality 
bin. The i? C oii uncertainties have been calculated 
by evaluating the changes in i? co n due to variations 
of the minimum-bias trigger efficiency, parameters 
of the Glauber calculation, and parameters in the 
modelling of the E£'^ Cal distribution [25 ]. The i? coU 
values and uncertainties are also reported in Ta- 
ble □ 



4. Monte Carlo samples 

Three Monte Carlo (MC) samples [25] were used 
for the analysis in this Letter. A total of 1 mil- 
lion simulated minimum-bias Pb+Pb events were 
produced using version 1.38b of the HIJING event 
generator [30]. HIJING was run with default pa- 
rameters except for the disabling of jet quenching. 
To simulate the effects of elliptic flow in Pb+Pb 
collisions, a parameterized centrality-, rj- and pt- 
dependent cos 2(f> modulation based on previous 
ATLAS measurements |25j was imposed on the par- 
ticles after generation [31]. The detector response 
to the resulting HIJING events was evaluated us- 
ing GEANT4 32J configured with geometry and 
digitization parameters matching those of the 2010 
Pb+Pb run. 

A "MC overlay" data set, intended specifically 
for evaluating jet performance, was obtained by 
overlaying GEANT4-simulated y/s^ = 2.76 TeV 
pp hard-scattering events on the HIJING events 
described above. The pp events were obtained 
from the ATLAS MC09 tune [33J of the PYTHIA 
event generator [31]. One million PYTHIA hard- 
scattering events were generated for each of five in- 
tervals of pt, the transverse momentum of outgoing 
partons in the 2 — > 2 hard-scattering, with bound- 
aries 17, 35, 70, 140, 280 and 560 GeV. The pp events 
for each j>t interval were overlaid on the same sam- 
ple of HIJING events. 

A smaller sample of "data overlay" events was 
produced by overlaying 150k GEANT4-simulated 
PYTHIA pp events onto 150k minimum-bias 
Pb+Pb data events recorded during the 2011 LHC 
Pb+Pb run. Due to the different detector condi- 
tions in the 2010 and 2011 runs, the data overlay 
events cannot provide the corrections required for 
this analysis. However, they provide a valuable test 
of the accuracy of HIJING in describing the under- 
lying event. 
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5. Jet reconstruction 



separately for each calorimeter layer according to 



Calorimeter jets were reconstructed from Arj x 
A(f) = 0.1 x 0.1 towers using the anti-fej algorithm 
|21) in four-vector recombination mode with anti-/cf 
distance parameters R — 0.2, 0.3, 0.4 and 0.5. The 
tower energies were obtained by summing energies, 
calibrated at the electromagnetic energy scale [55] . 
of all cells in all layers within the r\ and cj> boundaries 
of the towers. Cells that span tower boundaries had 
their energy apportioned by the fraction of the cell 
contained within a given tower. The jet measure- 
ments presented here were obtained by performing 
the anti-fct reconstruction on the towers prior to 
underlying event (UE) subtraction and then evalu- 
ating and subtracting the UE from each jet at the 
calorimeter cell level. The subtraction procedure 
calculates a per-event average UE energy density 
excluding contributions from jets and accounting 
for effects of elliptic flow modulation on the UE [36 . 
The UE estimation and subtraction was performed 
using a two-step procedure that was identical for all 
jet radii. 

A first estimate of the UE average transverse en- 
ergy density, Pi(^), was evaluated in 0.1 intervals of 
r\ from all cells in each calorimeter layer, i, within 
the given rj interval excluding those within "seed" 
jets. In the first subtraction step, the seeds are de- 
fined to be R = 0.2 jets containing at least one 
tower with Et > 3 GeV and having a ratio of 
maximum tower transverse energy to average tower 
transverse energy, EJf ax / (Et) > 4. Elliptic flow in 
Pb+Pb collisions can impose a 2i>2 cos [2(</> — ^2)] 
modulation on the UE. Here, V2 is the second coef- 
ficient in a Fourier decomposition of the azimuthal 
variation of the UE particle or energy density, and 
the event plane angle, ^2, determines the phase 
of the elliptic modulation. Standard techniques 
[25l [36] were used to measure VP 2 , 



/ w k E T k sin (20 fc )\ 



*2 = r tan~ 



^Wfc^TfeCOS (2(/) k ] 



(1) 



where k runs over cells in the FCal, (f>k represents 
the cell azimuthal angle, and Wk represent per-cell 
weights empirically determined to ensure a uniform 
$2 distribution. An ry-averaged V2 was measured 



v 2i 



^£ Tj cos [2(^-* 2 )] 



(2) 



where j runs over all cells in layer i. The UE- 
subtracted cell transverse energies were calculated 
according to 



sub 



= E Tj -Aj pifa) (1 + 2v 2l cos [2 (<k - * 2 )]) , 

. (3) 

where -Erjj rjj, <fij an d Aj represent the cell Et, rj 
and 4> positions, and area, respectively for cells, j, in 
layer i. The kinematics for R — 0.2 jets generated 
in this first subtraction step were calculated via a 
four-vector sum of all (assumed massless) cells con- 
tained within the jets using the Et values obtained 
from Eq. [3j 

The second subtraction step starts with the def- 
inition of a new set of seeds using a combination 
of R = 0.2 jets from the first subtraction step with 
Et > 25 GeV and track jets (defined below) with 
Pt > 10 GeV. Using this new set of seeds, a new 
estimate of the UE, p'^rf), was calculated exclud- 
ing cells within AR = 0.4 of the new seed jets, 
where AR = V'feeii - Vjct) 2 



t ) 2 . New 



'cell — Vict) 

V2i values, v^, were also calculated excluding all 
cells within A77 = 0.4 of any of the new seed jets. 
This exclusion largely eliminates distortions of the 
calorimeter V2 measurement in events containing 
high-px jets. The background subtraction was then 
applied to the original cell energies using Eq. [3] but 
with pi and V2i replaced by the new values, p'^rj) 
and V2i- New jet kinematics were obtained for all 
jet radii from a four-momentum sum of cells within 
the jets using the subtracted cell transverse ener- 
gies. Jets generated in this second subtraction step 
having Et > 20 GeV were recorded for subsequent 
analysis. 

A correction of typically a few per cent was ap- 
plied to the reconstructed jets to account for incom- 
plete exclusion of towers within jets from the UE es- 
timate due, for example, to differences in direction 
between the seeds and the final jets. This correc- 
tion was validated by applying the full heavy ion 
jet reconstruction procedure to 2.76 TeV pp data 
collected by ATLAS in March 2011. The recon- 
structed jets were compared, jet-by-jet, to those 
obtained from the pp jet reconstruction procedure. 
After this last correction for incomplete exclusion of 
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jets from the background, the energy scales of the 
heavy ion and pp reconstruction procedures agreed 
to better than 1% for E T > 25 GeV. A final cor- 
rection depending on the jet 77, Et, and R was ap- 
plied to obtain the correct hadronic energy scale for 
the reconstructed jets. The calibration constants 
were derived separately for the four jet radii using 
the same procedure applied to pp jet measurements 

m 

In addition to the calorimeter jet reconstruction, 
track jets were reconstructed using the anti-fc t al- 
gorithm with R = 0.4 from charged tracks that 
have a good match to the primary vertex and that 
have pt > 4 GeV. This threshold suppresses con- 
tributions of the UE to the track jet measurement. 
Specifically, an R = 0.4 track jet has an estimated 
likelihood of including an uncorrelated pt > 4 GeV 
charged track of less than 4% in the 0-10% central- 
ity bin. 

The fluctuating UE in Pb+Pb collisions can po- 
tentially produce reconstructed jets that do not 
originate from hard-scattering processes. In the re- 
mainder of this Letter such jets are referred to as 
"underlying event jets" or UE jets. A requirement 
that calorimeter jets match at least one track jet 
with pt > 7 GeV or an EM cluster reconstructed 
from cells in the electromagnetic calorimeter [57] 
with pt > 7 GeV was applied to exclude UE jets. 
The matching criterion for both track jets and EM 
clusters is that they lie within AR = 0.2 of the 
jet. Applying this matching requirement provides 
a factor of about 50 rejection against UE jets while 
inducing an additional px-dependent inefficiency in 
the jet measurement. To accommodate the use 
of track jets in the UE jet rejection, the jet mea- 
surements presented here have been restricted to 
\r)\ < 2.1. The total number of jets above pt thresh- 
olds of 40 GeV and 100 GeV in the data sample after 
event selection, UE jet rejection, and the |?y| < 2.1 
cut have been applied is shown in Table [2] for the 
most central and peripheral bins. 

6. Performance of the jet reconstruction 

The primary evaluation of the combined per- 
formance of the ATLAS detector and the analy- 
sis procedures described above in measuring un- 
quenched jets was obtained using the MC over- 
lay sample. In that MC sample, the kinematics of 
the reference PYTHIA generator-level jets (here- 
after called "truth jets") were reconstructed from 
PYTHIA final-state particles for R = 0.2,0.3,0.4 



Table 2: Total number of jets in the data set with px > 
40 GeV and p T > 100 GeV in the 0-10% and 60-80% cen- 
trality bins after all event selection criteria, UE jet rejection, 
and the \rj\ < 2.1 cut have been applied. 



R 


p T > 40 GeV 


p T > 100 GeV 


0-10% 


60-80% 


0-10% 


60-80% 


0.2 


112 333 


8068 


2308 


162 


0.3 


287 153 


12 629 


3534 


222 


0.4 


543 444 


15 964 


4974 


277 


0.5 


710 158 


18 573 


7586 


307 



and 0.5 using the same techniques as applied in pp 
analyses [35] . Separately, the presence and approxi- 
mate kinematics of HIJING-generated jets were ob- 
tained by running R — 0.4 anti-fc t reconstruction on 
final-state HIJING particles having pt > 4 GeV. 
Accidental overlap of jets from unrelated hard- 
scattering processes may occur at non-negligible 
rates in the data due to the geometric enhancement 
of hard-scattering rates in Pb+Pb collisions. How- 
ever, for the purposes of this Letter, the resulting 
combined jets are considered part of the physical jet 
spectrum and not a result of UE fluctuations. Then, 
to prevent the overlap of PYTHIA and HIJING jets 
from distorting the jet performance evaluated rela- 
tive to PYTHIA truth jets, all PYTHIA truth jets 
within AR = 0.8 of a p T > 10 GeV HIJING jet 
were excluded from the analysis. 

Following reconstruction of the overlaid MC 
events using the same algorithms that were ap- 
plied to the data, PYTHIA truth jets passing the 
HIJING-jet exclusion were matched to the clos- 
est reconstructed jet of the same R value within 
AR = 0.2. The resulting matched jets were used 
to evaluate the jet energy resolution (JER) and the 
jet energy scale (JES). The jet reconstruction effi- 
ciency was defined as the fraction of truth jets for 
which a matching reconstructed jet is found. The 
efficiency was evaluated both prior to (e) and fol- 
lowing (e') UE jet rejection. For all three perfor- 
mance measurements, the different pt MC overlay 
samples were combined using a weighting based on 
the PYTHIA cross-sections for each pt range. 

Figure [TJ shows a summary of the ATLAS Pb+Pb 
jet reconstruction performance for R = 0.2 and 
R = 0.4 jets in central (0-10%) and peripheral (60- 
80%) collisions. The (fractional) JER was char- 
acterized by cr[A£ , T ]/£; T ruth , where a[AE T ] is the 
standard deviation of the AE T = E r T cc - # T ruth 
distribution and where _E T CC and E^ are the re- 



5 




a> 0.8 : 
o 

3= 0.6 : 
LU 

0.4 : 



Uj 0.4 

' H 

LU 

^ 0.3 
O 

o 



- § 0.2 



• E, 0-1 0% 

e, 60-80% 

o e', 0-10% 

♦ £', 60-80% 



o 



5 



LU 



50 



100 150 200 250 



300 



Ef th [GeV] 



ATLAS simulation 

anti-/c, R = 0.4 

— e— o [A£J/£™ th + fit, 0-10% 
o [A£J/£" u,h + fit, 60-80% 
<A£ >/£" u,h , 0-10% 




• £,0-10% 

£, 60-80% 

o £',0-10% 

♦ £', 60-80% 



50 100 150 200 250 300 
Ef* [GeV] 



Figure 1: Results of MC evaluation of jet reconstruction performance in 0-10% and 60-80% collisions as a function of truth 
jet Et for R = 0.2 (left) and R = 0.4 (right) jets. Top: jet energy resolution a[AE T ]/E^ uth and jet energy scale closure, 
(AEt) I E^ ruth . Solid curves show parameterizations of the JER using Eq. 4 Bottom: Efficiencies, e and e' , for reconstructing 
jets before and after application of UE jet removal (see text for explanation), respectively. 



constructed and truth jet E^ values, respectively. 
The JES offset or "closure" was evaluated from the 
mean fractional energy shift, (AET)/E^ uth . 

The JER was found to be well described by a 
quadrature sum of three terms, 

a[AE T ] _ a b 

ptruth - / 7 = ® ptruth ® C ' W 
-C/fp / ptruth _CVrp 

where a and c represent the usual sampling and con- 
stant contributions to calorimeter resolution. The 
term containing b describes the contribution of un- 
derlying event fluctuations, which do not depend on 
jet E^. Results of fitting the Ef dependence of the 
JER according to Eq. |4j using methods described 
below, are shown with curves in Fig. [I] 

The jet reconstruction efficiency decreases with 
decreasing jet Et for Et < 50 GeV. The decrease 
starts at larger Et and decreases more rapidly for 
larger jet radii and in more central collisions. The 
inefficiency results primarily from the finite JER 
which causes jets with £^ ruth > 20 GeVto be mea- 
sured with i^ oc < 20 GeV. The UE jet rejection 
causes an additional loss of jets but in a manner 
that reduces the centrality dependence of the inef- 
ficiency. 



The accuracy of the MC overlay studies described 
above was evaluated using the data overlay sam- 
ple analyzed using the same procedures that were 
applied to the MC overlay sample. The analysis 
yielded results for the JER, JES, and efficiency con- 
sistent with the MC overlay sample, although the 
JER in the data overlay sample was found to be 
slightly better than in the MC overlay sample. The 
JES in the data overlay sample was found to agree 
between peripheral and central collisions to better 
than 1% for R = 0.4 jets, and the reconstruction ef- 
ficiency was found to differ by less than 5% on the 
rise of the efficiency curve. 

A data-driven check of the HIJING description of 
UE fluctuations was performed by evaluating dis- 
tributions of EM-scale summed E^ in rectangular 
groups of towers within the interval \rj\ < 2.8. The 
groups were chosen to match the areas of jets used 
in this analysis: 3x4 and 7 x 7 for R = 0.2 and 
R = 0.4 jets, respectively. No attempt was made to 
exclude jets from the fluctuation analysis. The dis- 
tributions of E T x4 and E T xT , the E£; T for 3 x 4 and 
7x7 groups of towers, are shown in Fig.[2]for a nar- 
row range of ££f Cal , 3.4 < E£:| Cal < 3.5 TeV, that 
lies within the 0-1% centrality interval. These dis- 
tributions have mean values, (£t X4 ) = 26 GeV and 
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Figure 2: Top: Representative distributions of Kj, x * - (Ej, } (left) and E^T - {E^ x '} (right) (see text for definitions) 
for data (points) and MC (filled histogram) for Pb+Pb collisions with 3.4 < SB^ Cal < 3.5 TeV. The vertical lines indicate 
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(E' x ') = 0. Bottom: Standard deviations of the Ei x4 and E' x ' distributions, <t[J5j.* ] and 



t\E\ 



respectively, in data and HIJING MC sample as a function of E_E£ Cal . 



(E 7 T x7 ) = 105 GeV, subtracted and, thus, m prin- 
ciple represent the distribution of the residual con- 
tributions of the UE to jet energies after subtrac- 
tion. However, the high tails of the distributions 
can be attributed to the presence of jets, which are 
not part of the UE. The corresponding distribu- 
tions obtained from the HIJING MC sample, but 
with (E^ xi ) and (E 7 ^ 7 ) obtained from data, are 
shown in Fig. [2] with filled histograms. 

The shapes of the MC and data distributions in 
Fig. [2] (top) are very similar, but the MC result 
slightly over-predicts the positive fluctuations for 
all collision centralities. In central collisions the 
MC result also slightly over-predicts the size of neg- 
ative fluctuations. In contrast, for non-central col- 
lisions (not shown here) the data has a broader 
distribution of negative fluctuations than the MC 
sample. These observations are demonstrated by 
Fig. [2] (bottom) which shows the standard devia- 
tions of the £^ x4 and E 7 T x7 distributions, a[E^ xA ] 
and (j[£t X7 ], as a function of E£^ Cal , obtained 
from both the data and the MC sample. The data 
and MC distributions have similar trends, but the 
MC cr[£^ x4 ] and u[i^ x7 ] values are larger in cen- 
tral collisions by 2.5% and 5%, respectively. In non- 



central collisions, the broader spectrum of negative 
fluctuations in the data causes af-E^ 4 ] and alE 7 ^ 7 ] 
to exceed the corresponding quantity in the HIJING 
MC sample by approximately the same percentages. 

Consistency between the results of the fluctua- 
tion analysis and the evaluation of the JER de- 
scribed above has been established by fitting the 
Et dependence of the JER with the functional form 
given by Eq.[4]with fixed b values obtained from the 
fluctuation analysis. The b values for a given jet ra- 
dius were determined by taking the standard devi- 
ation of the Si?T distribution for the corresponding 
tower group averaged over centrality and corrected 
to the hadronic energy scale. The resulting b val- 
ues for R = 0.2(0.4) jets are 5.62(12.45) GeV and 
1.15(2.58) GeV for the 0-10% and 60-80% central- 
ity bins respectively. The parameters a and c ob- 
tained from the fits are found to be independent 
of centrality within fit uncertainties, as expected, 
and to have values a = 1.0(0.8), c = 0.07(0.06) for 
R = 0.2(0.4) jets with E T expressed in GeV. The 
accuracy of the fits in describing the Et dependence 
of the JER is demonstrated by the curves showing 
results for R = 0.2 and R = 0.4 jets in Fig. [lj 

The contribution of UE jets to the measured 
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jet spectrum after UE jet rejection is estimated 
to be < 3% approximately independent of jet pt 
for 40 < pt < 60 GeV and less than 1% for 
Pt > 60 GeV. This estimate was obtained by eval- 
uating the rate of reconstructed jets in the HIJING 
MC sample which were not matched to HIJING 
truth jets and correcting for missing truth jets due 
to the pt > 4 GeV requirement applied in the HI- 
JING truth jet reconstruction. 

7. Jet spectra and unfolding 

Though jet reconstruction performance is nat- 
urally evaluated in terms of jet Et, the physics 
measurements in this Letter were performed as a 
function of pt directly calculated from the jet four- 
momentum. The typical masses of the jets are suf- 
ficiently small that Et ~ Pt holds over the range 
of measured pt for all jet radii. The measured pt 
spectra of reconstructed jets passing UE jet rejec- 
tion and having \rj\ < 2.1 were evaluated for each 
centrality bin using logarithmic pt bins spanning 
the range 38 < pt < 210 GeV. The correlations 
within and between pt bins arising from multi-jet 
events were quantified by the covariance, Cy, be- 
tween the number of jets measured in two bins, i 
and j. The measured i?cp was calculated as 



/ N^(pt) \ 



pmcas 



(Pt) 



Dccnt 
^coll 



1 v evt 



iVfet 



V 



]\r60 80 
1 v evt 
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where iV-^. 11 * represents the measured jet yield in a 
given pt and centrality bin, and iV™"* and N®®^~ 80 
are the number of Pb+Pb collisions within the cho- 
sen and peripheral reference centrality intervals, re- 
spectively. Results for i?™p as |o-io obtained from the 
measured spectra are shown in Fig. [3] for R = 0.2 
and R = 0.4 jets. The i?^ as | -io for R = 0.2 jets 
is approximately equal to 0.5 over the measured 
p T range. The flg| as |o-io for R = 0.4 and R = 0.2 
jets are consistent for pt > 120 GeV, but at lower 
Pt, the R — 0.4 i?cp as |o-io increases relative to the 
R = 0.2 values. The difference between R = 0.2 and 
R = 0.4 -Rcp as |o-io values can be mostly attributed 
to the difference in the size of the UE fluctuations 
for R = 0.2 and R = 0.4 jets shown in Fig. [I] The 
larger JER for R = 0.4 jets produces greater up- 
ward migration on the steeply falling jet pt spec- 
trum in central collisions than in peripheral col- 
lisions, thus enhancing the measured Rcp- The 
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Figure 3: Top: Measured and corrected i?cp values for 
the 0—10% centrality bin as a function of jet px for R = 0.4 
and _R = 0.2 jets. Bottom: Ratio of corrected to measured 
i?Cp values for both jet radii. The error bars on the points 
represent statistical uncertainties only. 



drop in the R = 0.4 i?™p as |o-in at low p T is due 
to the decrease in jet reconstruction efficiency be- 
tween 60-80% and 0-10% centrality bins which, as 
noted above, largely results from the worse JER in 
central collisions. 

To remove the effects of the bin migration, the jet 
spectra were unfolded using the singular value de- 
composition (SVD) technique [38] as implemented 
in RooUnfold The MC overlay samples were 

used to populate a response matrix, A, which de- 
scribes the transformation of the true jet spec- 
trum, x, to the observed spectrum, b, according 
to b — Ax. The truth and reconstructed jet pt 
were obtained from the MC overlay sample using 
the methods described in Section [6j and[5j respec- 
tively, and the selection and matching of truth and 
reconstructed jet pairs was performed as described 
in Section|6] Using the weighting method suggested 
in Ref. [38], the unfolded spectrum is expressed as 
a set of weights w multiplying the input spectrum 
(xi n i) used to produce A. The SVD method ex- 
presses the solution for w in terms of a least-square 
minimization problem that includes a prescription 
for regularizing the amplification of statistical fluc- 
tuations of the data that would result from the di- 
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rect inversion of A. The regularization is controlled 
by a parameter r such that contributions from sin- 
gular values Sk of the unfolding matrix with Sk < T 
are suppressed. Inclusion of the pr-dependent re- 
construction efficiency in the response was found to 
strongly affect the spectrum of singular values of 
the matrix defining the SVD problem, so the effi- 
ciency correction was applied separately following 
the unfolding. The spectrum of MC truth jets was 
reweighted to provide a smooth, power-law initial 
spectrum, x- m i oc e'(pT)/p T , wnere the power index 
was chosen to be n = 5. An analysis of the optimal 
regularization in the SVD unfolding following the 
methods of Ref. [38] indicated that a regularization 
parameter fixed by the fifth singular value (r = s|) 
of the SVD matrix was appropriate for all central- 
izes and all R values. The statistical uncertainties 
in the SVD unfolding due to statistical errors on the 
input spectrum were evaluated using the pseudo- 
experiment technique with 1000 separate stochastic 
variations of the input spectrum based on the full 
covariance matrix. The contributions of statistical 
fluctuations in the response matrix, A, were simi- 
larly evaluated using an equal number of stochastic 
variations of the response matrix. The two contri- 
butions to the statistical uncertainty were combined 
in quadrature. 

Potential biases in the unfolding procedure were 
evaluated using two different methods. Each un- 
folded spectrum was re-folded with its correspond- 
ing response matrix and compared to the measured 
spectrum for self-consistency. In general, regular- 
ization can introduce differences between re-folded 
and measured spectra on the scale of statistical 
uncertainties on the measured spectra, while over- 
regularization can produce larger, systematic differ- 
ences. For all of the unfolded spectra, the re-folding 
procedure yielded a typical difference between mea- 
sured and re-folded spectra comparable to the sta- 
tistical uncertainties on the measured spectra. A 
separate check was performed by unfolding the re- 
constructed MC spectrum for each centrality bin 
and each jet radius and comparing to the original 
MC truth jet spectrum. For this purpose, the MC 
data sets were divided in half and reconstructed 
spectra and response matrices were generated sep- 
arately from each set. The unfolded and truth 
MC jet spectra typically agreed to better than 2%, 
though for the 0-10% centrality bin and for R = 0.4 
and 0.5 jets, differences as large as 5% were ob- 
served in the lowest pt bins. These differences are 
covered by the unfolding systematic uncertainties 



described below. 

The corrected Rcp was evaluated according to 
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where iVjet represents the unfolded number of jets 

,i i . i /cent j /60— 80 ,i 

m the pt bin, and e and e are tne px _ 

dependent jet reconstruction efficiencies after UE 
jet rejection for the indicated centrality bins. Fig- 
ure [3] shows the comparison of the corrected and 
measured Rcp values as a function of jet pt for 
R = 0.2 and R = 0.4 jets in the 0-10% centrality 
bin. The unfolding has little effect on the R = 0.2 
Rcp due to the good energy resolution (relative to 
larger radii) for R = 0.2 jets even in central col- 
lisions. For the R = 0.4 jets, Rqp is reduced by a 
factor of about two at the lowest pr values included 
in the analysis and is only slightly modified at the 
highest pt- Because the unfolding provides a non- 
local mapping of the input jet pr spectrum onto 
the unfolded spectrum, the statistical uncertainties 
in the unfolded spectra have significant correlations 
between bins, and there is not a direct relationship 
between the statistical errors in the input spectrum 
and the unfolded spectrum. The regularization of 
the unfolding also suppresses statistical fluctuations 
in the unfolded spectrum, but the statistical uncer- 
tainties in the measured spectrum also contributes 
to the systematic uncertainties from the unfolding 
procedure. 

8. Systematic uncertainties 

Systematic uncertainties in the Rqp measure- 
ment can arise due to errors on the jet energy scale 
(JES), the jet energy resolution (JER), jet finding 
efficiency, the unfolding procedure, and the i? co n 
values. Uncertainties in jet E^ and px are assumed 
to be equal (i.e. Sp^ = SEt). Uncertainties in the 
JES and the JER influence the unfolding of the jet 
spectra. The resulting systematic uncertainties on 
the i?cp values (5R^p) were evaluated by produc- 
ing new response matrices according to the proce- 
dures described below, generating unfolded spectra 
from these matrices, and calculating new Rcp val- 
ues. The resulting changes in the Rcp values were 
taken to be estimates of SR^p. For uncertainties 
fully correlated in centrality, SR^p was evaluated 
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by simultaneously varying the chosen centrality bin 
and the 60-80% bin, while for other uncertainties, 
the chosen centrality bin and 60-80% centrality 
bins were varied separately and the variations in 
Rcp combined in quadrature. 

Overall JES uncertainties common to the differ- 
ent centrality bins cancel in the ratio of the spectra 
in Rcp, but ccntrality-dependent JES errors will 
produce a systematic shift in i?cp- Studies using 
the MC overlay sample discussed in Section [6] in- 
dicate a maximum difference in JES between the 
0-10% and 60-80% centrality bins for the jet pt 
range included in this analysis of 0.5%, 1%, 1.5% 
and 2.5% for R = 0.2,0.3,0.4 and 0.5, respectively. 
Studies were also performed with the data overlay 
sample using an identical procedure as that applied 
to the MC overlay sample. The JES evaluated in 
the data overlay was found to agree between the 
0-10% and 60-80% centrality bins to better than 
1%, which is better than the agreement found in 
the MC overlay sample. 

Independent evaluations of a possible central- 
ity dependence of the calorimeter JES were per- 
formed by matching track and calorimeter jets in 
both the data and the MC overlay sample. The 
track jets provide a common reference for evalu- 
ating calorimeter jet response that is insensitive 
to the UE. The average calorimeter jet E^ was 
evaluated as a function of matching track jet pr, 
(E^ d °)(p t ^ ct ), for different centrality bins. In the 
data, for p£ kjet > 50 GeV, the (£^ al °) values were 
found to be consistent across all centrality bins to 
better than 3%. Accounting for a slight centrality 
dependence seen in the MC overlay sample, the 0- 

< 



trkjet 



10% and 60-80% bins agree to 2%. For p™ 
50 GeV, R- and centrality-dependent differences of 
up to 4% (for R — 0.5) are observed between data 
and MC overlay results for (E^°)(p^ kiet ). This 
study provides a stringent constraint on changes 
in calorimeter response for jets affected by quench- 
ing and justifies the use of unquenched jets from 
PYTHIA in evaluating the jet performance and re- 
sponse matrices. 

Based on the combination of the studies de- 
scribed above, the systematic uncertainties on the 
centrality dependence of the JES for the 0-10% cen- 
trality bin and for calorimeter jet px > 70 GeV were 
estimated to be 0.5%, 1%, 1.5% and 2.5%, respec- 
tively, for R = 0.2,0.3,0.4 and 0.5 jets. At lower 
Pt, the assigned systematic uncertainties increase 
linearly with decreasing px such that they double 
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Figure 4: Contributions to the relative systematic uncer- 
tainty on the Rqp from various sources for the R = 0.4 
anti-fct jets in the 0-10% centrality bin. The k ± 1 curves 
denote the uncertainty due to the choice of regulariza- 
tion parameter obtained by unfolding with the fourth and 
sixth singular values. A constant 5% systematic uncer- 
tainty on the jet reconstruction efficiency is assigned for 
Pt < 100 GeVonly. The 11% uncertainty in the determi- 
nation of i? C oii is indicated with a shaded box and is px- 
independent. 



in size between 70 GeV and 38 GeV. For other cen- 
trality bins, the systematic errors on the central- 
ity dependence of the JES decrease smoothly from 
central to peripheral collisions. The resulting 5R^p 
values were evaluated using new response matrices 
generated by scaling the reconstructed px to ac- 
count for the above-quoted JES uncertainties. The 
JES systematic uncertainty is assumed to be fully 
correlated between different centrality bins and dif- 
ferent R values. 

Systematic uncertainties in the JER due to in- 
accuracies in the MC description of the UE fluctu- 
ations were evaluated using results of the fluctua- 
tion analysis described above. The effects of those 
inaccuracies were evaluated by rescaling the per-jet 
Apt = _Px c — p^ uth values obtained from the MC 
study by factors that cover the differences between 
data and MC result. For each centrality and jet 
radius, a modified value of the b parameter in Eq. [4] 
was evaluated and used to obtain new JER values, 
<j'[AEt] from Eq. |4j Then a rescaled Apx was ob- 
tained from 



Apx' = Apx 



(7) 



Since the discrepancies between the MC and the 
data were observed to be different for positive and 
negative fluctuations, the rescaling was applied sep- 
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arately for positive and negative Apx. 

The YjEt values in the MC study were found 
to have larger positive fluctuations than those in 
the data for all centralities by approximately 2.5%, 
2.5%, 5%, and 7.5% for R = 0.2,0.3,0.4 and 0.5 
jets, respectively, so for positive Apx b was reduced 
by these percentages. For the 0-10% centrality bin, 
the negative fluctuations were also larger in the MC 
study than in the data by the same approximate 
percentages, so for central collisions the same, mod- 
ified b value was used for negative Apx- For all 
other centrality bins, the negative fluctuations in 
the data were larger than in the MC by approxi- 
mately twice the above-quoted percentages. Thus, 
for those centralities, the modified b values were ob- 
tained for negative Apx by increasing b by 5%, 5%, 
10%, and 15%, respectively, for R = 0.2,0.3,0.4 
and 0.5 jets. 

New response matrices were generated using the 
calculated Apx' values according to p™ c ' = p^ uth + 
Apx', and these modified response matrices were 
used to estimate the JER systematic uncertain- 
ties following the procedure described above. The 
systematic uncertainty on the spectra due to the 
JER for the 0-10% centrality bin was taken to be 
one-sided as all evaluations indicate that the MC 
simulations slightly overestimate UE fluctuations. 
Asymmetric errors were obtained for the other cen- 
trality bins by applying the positive and negative 
AEt scalings separately. The JER systematic un- 
certainties were assumed to be fully correlated be- 
tween different jet R values but uncorrelated be- 
tween different collision centralities, so the uncer- 
tainties on the spectra were combined in quadrature 
in evaluating SR^p. The conservative assumption 
that the JER uncertainties are fully uncorrelated 
between different centrality bins is based on the ob- 
servation that the differences between data and the 
HIJING MC sample in the fluctuation analysis are 
not the same for all centralities. 

The systematic uncertainties associated with the 
non-UE contributions to the JER (described by the 
a and c terms in Eq. [4]) were evaluated following 
procedures used by ATLAS in previous pp jet mea- 
surements }40j . New response matrices were gener- 
ated by applying an additional stochastic smearing 
to the Apx values, and the systematic uncertainty 
was obtained by applying the procedure described 
above. 

Systematic uncertainties on Rcp due to the un- 
folding were evaluated by changing the power in- 
dex (n) in the functional form for a;i n ; by ±0.5 and 



by varying the regularization parameter. The ±0.5 
change in the power law index was chosen because 
it produces a spectrum that changes relative to the 
default Xi n i over the measured px range by a fac- 
tor of about two - the typical suppression observed 
in central collisions. Thus, it covers the possibility 
that the true i?cp could increase to one or decrease 
to 0.25 over the measured px range. To evaluate 
the potential systematic uncertainty due to regu- 
larization, the unfolding was performed with regu- 
larization parameters obtained from the fourth and 
sixth singular values of the unfolding matrix, r = s\ 
and r = s§. Systematic uncertainties on the spec- 
tra were determined from the differences in the un- 
folded spectra. The resulting SR s ^p values were ob- 
tained assuming that the regularization uncertain- 
ties on the two spectra are uncorrelated. 

The systematic uncertainty on the efficiency cor- 
rection was evaluated by comparing MC overlay 
and data overlay samples where differences less 
than 5% were observed on the "turn on" part of 
the efficiency curve. A 5% uncertainty due to 
the efficiency correction was applied to Rcp for 
Pt < 100 GeV in the four most central bins. 

The different contributions to the total SRqp are 
shown in Fig. [I] for R = 0.4 jets in the 0-10% cen- 
trality bin. The JES and Xi n i uncertainties are ap- 
proximately independent of px, while the JER un- 
certainty decreases with increasing px. The regu- 
larization uncertainty grows with increasing px due 
to the poorer statistical precision of the high-px 
points. The systematic uncertainties for the other 
radii show similar px and centrality dependence, 
with the JES and JER uncertainties increasing with 
jet radius as expected. 

9. Results 

Figure [5] shows the Rqp values obtained for 
R = 0.2 and R = 0.4 jets as a function of px in 
four bins of collision centrality with three different 
error contributions: statistical uncertainties, par- 
tially correlated systematic uncertainties, and fully 
correlated uncertainties. The Rqp values for all 
centralities and for both jet radii are observed to 
have at most a weak variation with px. For the 
0-10% centrality bin the Rqp values for both jet 
radii show a factor of about two suppression in the 
1/ A co n-scaled jet yield. For more peripheral colli- 
sions, Rqp increases at all jet px relative to central 
collisions, with the Rqp values reaching 0.9 for the 
50-60% centrality bin. A more detailed evaluation 
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Figure 5: Rcp values as a function of jet px for R = 0.2 (left) and R = 0.4 (right) anti-fct jets in four bins of collision centrality. 
The error bars indicate statistical errors from the unfolding, the shaded boxes indicate unfolding regularization systematic errors 
that are partially correlated between points. The solid lines indicate systematic errors that are fully correlated between all 
points. The horizontal width of the systematic error band is chosen for presentation purposes only. Dotted lines indicate 
Rqp = 0.5, and the dashed lines on the top panels indicate i?cp = 1- 
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presentation purposes only. Dotted lines indicate Rcp = 0.5, and the dashed lines on the top panels indicate Rcp = 1. 



of the centrality dependence of Rcp for R = 0.4 jets 
is presented in Fig. [6j which shows Rcp vs A/p art for 
six jet pt bins. Rqp decreases monotonically with 
increasing N part for all px bins. The lower p^ bins, 
for which the data are more statistically precise, 
show a variation of Rqp with N pai t that is most 
rapid at low A^ part . Trends similar to those shown 
in Figs. [5] and [6] are observed for all jet radii. 

The dependence of Rqp on jet radius is shown in 
Fig.[7]for the 0-10% centrality bin in four jet pt in- 
tervals (left) and for different centrality bins in the 
89 < pt < 103 GeV bin (right). For this figure, 
the shaded boxes indicate the combined contribu- 
tion of systematic uncertainties due to regulariza- 
tion, Xini, and efficiency, which are only partially 
correlated between points. All other systematic er- 
rors are fully correlated and are indicated by solid 
lines. The results in Fig. [7] show a weak variation 
of i?cp with R, that is nonetheless significant when 
taking into account the correlations in the errors 
between the different R values. 

To demonstrate this conclusion more clearly, 
Fig. [8] shows the ratio of Rqp values between R = 
0.3, 0.4 and 0.5 jets and R = 0.2 jets, R% P /Rgp, as 
a function of pt for the 0-10% centrality bin. When 
evaluating the ratio, there is significant cancella- 
tion between the correlated systematic uncertain- 
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Figure 8: Ratios of Rcp values between R = 0.3, 0.4 and 
0.5 jets and R, = 0.2 jets as a function of px m the 0-10% 
centrality bin. The error bars show statistical uncertainties 
(see text). The shaded boxes indicate partially correlated 
systematic errors. The lines indicate systematic errors that 
arc fully correlated between different px bins. 



ties. Statistical correlations between the jet yields 
for the different radii were evaluated in the mea- 
sured spectra and tracked through the unfolding 
procedure separately for the 0-10% and 60-80% cen- 
trality bins. Those correlations were then included 
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when evaluating the statistical errors on Rjip/Rjip 
shown in Fig. [§] The results in that figure indi- 
cate a significant dependence of Rcp on jet radius. 
For p T < 100 GeV the Rgp/Rgp 2 values for both 
R = 0.4 and R = 0.5 differ from one beyond the sta- 
tistical and systematic uncertainties. The deviation 
persists for R = 0.5 above 100 GeV. A similar, but 
weaker dependence is observed in the 10-20% cen- 
trality bin. In more peripheral bins, no significant 
radial dependence is observed. The differences be- 
tween i?cp values for the different jet radii increase 
with decreasing px, except for the lowest two px 
bins. However, direct comparisons of Rcp between 
different jet radii at low px should be treated with 
care as the same jets measured using smaller radii 
will tend to appear in lower px bins than when mea- 
sured with a larger radius. 

10. Conclusions 

This Letter presents results of measurements of 
the centrality dependence of jet suppression, char- 
acterized by the inclusive jet central-to-peripheral 
ratio, Rcp, in Pb+Pb collisions at 2.76 TeV per 
nucleon at the LHC. The measurements were per- 
formed over the px range 38 < px < 210 GeV for 
anti-fct jets of radii R = 0.2,0.3,0.4 and 0.5. The 
inclusive jet yield is observed to be suppressed by 
a factor of about two in central collisions relative 
to peripheral collisions with at most a weak px 
dependence to the suppression. The suppression 
varies monotonically with collision centrality over 
the measured px range and for all jet radii. The 
Rcp at fixed px is observed to vary with jet ra- 
dius increasing gradually from R — 0.2 to R = 0.5. 
That variation is most significant for px < 100 GeV 
where more than a 30% variation is observed. These 
results provide the first direct measurement of in- 
clusive jet suppression in heavy ion collisions. The 
substantial suppression of the jet yield observed 
at all px values complements the previous mea- 
surements of dijet transverse energy imbalance in 
Pb+Pb collisions at the LHC [IMS]. 
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